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ABSTRACT 

We present a simple dimensional argument to illustrate the impact of nonthermal support from turbulent velocity 
dispersion on the shape of the prestellar core mass function (CMF), precursor of the stellar initial mass function (IMF). 
The argument demonstrates the need to invoke such support to recover the Salpeter slope in the high-mass part of the 
CMF/IMF, whereas pure thermal support leads to a much steeper slope. This simple dimensional argument clearly 
highlights the results obtained in the complete Hennebelle-Chabrier theory of the IMF. 
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1. Introduction 

Understanding the shape of the stellar initial mass func- 
tion (IMF) has been an ongoing domain of research in 
astrophysics since the pioneering work of Salpeter (1955). 
Indeed, the IMF is the fundamental parameter governing 
stellar and galaxy evolution, and it is rooted in the very 
process of star formation. In the modern context of star 
formation, it is now widely accepted that star formation is 
a dynamical process, where turbulence plays a major role 
(see e.g Brunt 2003, MacLow & Klessen 2004, McKee & 
Ostriker 2007). Turbulence is injected on large scales in 
the ISM by various processes of stellar or galactic origin 
so that turbulent motions determine the structure of all 
density regimes in the interstellar gas and in star-forming 
molecular clouds (Elmegreen & Scalo 2004). In weak den- 
sity media, like molecular clouds, such large-scale super- 
sonic turbulence cascades to small scales by shocks. 

Turbulence is by nature a nonlinear process, developing 
from the strong distortions of the fluid velocity field due to 
advection (mathematically formalized by the (u ■ V) oper- 
ator in the hydrodynamic equations) . Interactions between 
these strong random distortions of the velocity field glob- 
ally act as a nonthermal source of support, stabilizing oth- 
erwise unstable scales. Often, and incorrectly for large-scale 
turbulence, such a support is denominated for convenience 
nonthermal " pressure" . 

The statistics of the gas density fluctuations created by 
turbulence is captured by the so-called probability density 
function (PDF). It is generally admitted that, in good ap- 
proximation, the PDF of turbulence resembles a lognor- 
mal form, as expected from the application of the central 
limit theorem to a hierarchical density field generated by 
multiplicative processes, such as shocks (Vazquez-Semadeni 
1994). Such a lognormal behaviour for density fluctuations 
has received observational support. Indeed, observationally, 
the 3D lognormal shape of the PDF is reflected by the pro- 



jected (2D) power spectrum column density of molecular 
clouds, measured from dust extinction maps (Kainulainen 
et al. 2009, Lombardi et al. 2010, Brunt et al. 2010). 
Determining the velocity field fluctuations is more prob- 
lematic. As one can normally only measure the velocity 
dispersion along the line of sight (cr z ), it is only possible 
to constrain the tangential velocity dispersion by assuming 
isotropy. 

In the case of compressible turbulence, the standard de- 
viation of the logarithm of the density fluctuations log(p/p), 
i.e. the width of the PDF, determined from numerical sim- 
ulations, is found to obey reasonably well the relation (e.g. 
Padoan & Nordlund 2002, Federrath et al. 2010) 

a 1 2 og(p) ^ln[l + (6Ai) 2 ], (1) 

where Ai is the hydrodynamical rms Mach number. This 
corresponds to a standard deviation o~ p ~ bM. for the den- 
sity fluctuation p/p. The factor b is found to be b as 0.5 
in 3D simulations for a mixture of solenoidal and com- 
pressive turbulence driving modes (Federrath et al. 2010, 
Schmidt et al. 2010), in agreement with observational deter- 
minations (Padoan, Jones & Nordlund 1997, Brunt 2010). 
Interestingly enough, the density PDF and relation ([T]) 
seem to remain weakly affected, for A4<10, in the case 
of supersonic MHD turbulence (Lemaster & Stone 2008, 
Kritsuk et al. 2011, Price et al. 2011). 

The first comprehensive theory of the IMF in the con- 
text of turbulent fragmentation of a molecular cloud has 
been formalized by Padoan & Nordlund (2002, PN). In 
this theory small-scale cascade of large-scale turbulence by 
shocks generates overdense regions that may become grav- 
itationally unstable, giving birth to the prestellar cores. 
A different theory has been proposed more recently by 
Hennebelle & Chabrier (2008, 2009, HC; see Hennebelle & 
Chabrier 2011 for a brief summary of these two theories). 
Although both theories predict that turbulence promotes 
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small-scale overdensities, triggering the formation of low- 
mass objects, in particular brown dwarfs, the two theories 
differ on several points. A fundamental difference, in par- 
ticular, concerns the role of nonthermal support on the de- 
termination of the IMF. On one hand, PN do not include 
turbulent support in their theory, and must invoke mag- 
netic field and MHD shock conditions to obtain the proper 
Salpeter slope for the IMF. In contrast, in the HC theory, 
random turbulent motions inside bound massive overdense 
regions, identified at the very early stages of star formation, 
counteract the action of gravity and thus stabilize these 
initial mass reservoirs - precursors of the gravitationally 
bound prestellar " cores" - which would otherwise have col- 
lapsed before collecting such an amount of mass. As such, 
these turbulent motions are acting like a source of non- 
thermal support (identified in the HC theory by the quan- 
tity the Mach number at the Jeans scale, see Eqs.(41) 
and (43) of HC08 ) against gravity. As demonstrated in 
HC08, this turbulent support is mandatory to obtain the 
proper Salpeter slope. However, by the time the (massive) 
cores/stars form out of the (massive) reservoirs, turbulence 
will have had ample time to dissipate. We come back to 
this important concept of turbulent "support" at the end 
of the paper. 

Conversely, a fundamental prediction of the HC theory 
is that, in the absence of non-thermal support, the slope of 
the IMF is significantly steeper than the Salpeter one. This 
result, hence the importance of the aforementioned source 
of nonthermal support in order to get the proper Salpeter 
slope for the IMF, has received some support from recent 
simulations aimed at exploring this issue (Schmidt et al. 
2010), even though more work is definitely needed to nail 
down this issue. 

In this short paper, we present a simple dimensional 
argument, based on the well-known scaling properties of 
turbulence, to demonstrate that indeed, turbulent support 
is needed to recover the proper Salpeter high-mass tail of 
the CMF/IMF. 

2. Scaling properties of turbulence 

The average statistical properties of a turbulent density and 
velocity field can be described by the scale dependence of 
the corresponding 3D power spectra in Fourier space, as 
described by power-law relationships, Psd(x) = P x (k) oc 
fc~™ x , where x denotes either the velocity (x = a), den- 
sity [x = p/p)) or logarithm of density (x = log(p/p)) 
fluctuations. The energy spectrum, E(k), which is the ki- 
netic energy per unit mass per unit wavenumber, is equal 
to the angular average of the power spectrum: E(k) = 
4Trk 2 {P 3 rj(k)) oc fc(~ rl+2 ' in 3D. It corresponds to the spec- 
trum of velocity fluctuations over scales smaller than a size 
i?, a % = J™ /R E(k)dk oc J^ /R k^ n+2 ^dk oc R n ~ 3 , where 

a\ describes the variance of the 3D velocity field, a(x, y, z), 
within a volume R 3 . In real space, the energy spectrum thus 
corresponds to the second-order rms velocity fluctuations 
measured on scale R, {cr\). 

The two well known limits of incompressible and pres- 
sureless turbulence correspond to the Kolmogorov and 
Burgers limits, with respective 3D velocity power spectrum 
indexes n = 11/3 and n — 4. Numerical simulations of 
compressible, shock-dominated turbulence under Mach val- 
ues typical of molecular cloud conditions tend to show that 



the power spectrum of turbulent velocity and logarithm of 
density fluctuations exhibit a power-law behaviour, with a 
similar index, namely n a ~ n iog(p) ~ 3.8 (Kritsuk et al. 
2007). 

3. Analytical dimensional argument 

According to the Hennebelle-Chabrier theory, in the ran- 
dom field of density fluctuations characteristic of the cloud's 
structure, initial pieces of fluid out of which prestellar cores 
and individual stars (or multiple star systems such as bi- 
naries) will eventually form, are identified as overdense re- 
gions that isolate themselves from the surrounding medium 
and start to contract under the action of gravity. The con- 
dition to select such a given marginally unstable piece of 
fluid of mass M and size i?, with internal velocity dis- 
persion aji = a, simply stems from the virial condition: 
2£*m + Ep 0t = 0, where Ekin ~ (l/2)pc 2 denotes the ki- 
netic energy and E pot ~ pC-^- the gravitational energy. 
Using dimensional analysis, this yields the relation 



Mo 2 ~G^ i.e. M~^. (2) 

In the limit where the velocity dispersion is dominated by 
thermal motions, and assuming nearly isothermal condi- 
tions, a does not depend on R and condition ([2]) yields 

M th oc R (THERMAL SUPPORT). (3) 

Conversely, in the limit where the velocity dispersion is 

(n-3) _. 

dominated by turbulent motions, a oc R 2 (see g2J), 
yielding for the marginal instability condition 

Mt U rb oc R {n - 2) oc R 2r > +1 (TURBULENT SUPPORT). (4) 

with r) — (n — 3 ) / 2 . The thermal case is recovered for n = 3 , 
i.e. 77 = 0. For n > 3, i.e. rj > 0, we note the stronger de- 
pendence of the mass M upon the size R in the presence 
of turbulent support, compared with the pure thermal case 
(Eq.Q). The complete relation is obtained in HC08 (their 
eqn.(26)) and illustrates that a bound stable region can 
achieve larger masses for a given size thanks to the non- 
thermal support at the characteristic Jeans length, natu- 
rally leading to a larger number of high mass reservoirs 
(then of massive prestellar cores), which would otherwise 
have collapsed long before. This impact of nonthermal sup- 
port on the mass-size relation of massive bound cores has 
been confirmed by recent numerical simulations (Schmidt 
et al. 2010, Fig. 7). 

Assuming a large enough multi-scale structured space of 
size L and dimension D, thus volume L D , the probability of 
having N gravitationally bound structures on a scale larger 
than R is given by the probability law V(Nr) = N(> R) oc 
R~ D . In Fourier space, this simply means that the num- 
ber of fluctuations of wavenumber k ~ 1/R is proportional 
to the volume of a fluctuation, i.e. dN(k) oc dk = kP~ l dk. 
The density probability, i.e. the probability to have a struc- 
ture of size e [R, R + dR], thus reads as 

^MocR-^. (5) 
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A homogeneous volume in 3 dimensions obviously cor- 
responds to I? = 3, i.e. a uniform density probability 
dN(k)/dk — constant, as assumed e.g. in Padoan & 
Nordlund (2002). 

From Eqs. Q and (JSJ) , the probability of having a bound 
structure of mass e [M (R) , M +dM (R+dR)} , which defines 
the CMF/IMF, reads 



dN 
dM 



oc M~ 



with a = 



1 + (D - 3) 2n + 4 + {D - 3) 



2^ + 1 



(6) 



As seen from Eq.©, for D = 3, pure thermal support 
(n = 3,7? = 0) yields a power-law IMF oc M~ 4 , 
substantially steeper than the Salpeter value, a — 2.3ifl 
In the two limiting cases of incompressible (Kolmogorov) 
(n = 11/3, T) = 1/3) and pressureless (Burgers) (n = 4, 77 = 
1/2) turbulence, the turbulent support yields a — 2.8 and 
a = 2.5, respectively. In molecular clouds, turbulence is as- 
sumed to scale according to the observed Larson (1981) ve- 
locity dispersion - size relation with a typical value 77 ~ 0.4- 
0.5. This value is well recovered for the aforementioned 
power spectrum index inferred from numerical simulations 
of supersonic turbulence, n = 3.8. According to eqn.©, 
this yields a slope for the CMF/IMF a — 2.66, between the 
Kolmogorov and Burgers values. 

There are observational suggestions, however, that 
molecular clouds have a fractal dimension (Elmegreen & 
Falgarone 1996), with D w 2.5-2.7 (Sanchez et al. 2005). 
Herschel observations of gravitationally bound prestellar 
cores, on the other hand, suggest a mass-size relation 
M oc Rl 3 , with /3 « 1-2 (Fig. 4 of Andre et al. 2010 and 
Konyves et al. 2010). Combined with a Salpeter mass dis- 
tribution, this yields dN/dR oc R~ {1+D) with D « 1.4- 
2.7. Although these results should be taken with caution, 
they suggest a structured space with D < 3 all the way 
from clouds to cores. This is consistent with turbulence 
being the dominant mechanism that initially determines 
the distribution of both unbound and bound structures. 
It arises, however, from two different properties of turbu- 
lence. As demonstrated in HC08, the structure of the un- 
bound clouds/clumps arises from the (log) density structure 
of turbulence, characterized by a lognormal PDF and a 
power spectrum index n' . On the other hand, the distri- 
bution of gravitationally bound structures is set up by the 
virial condition and thus involves the properties of the ve- 
locity structure of turbulence (see Eq. $Z§), characterized 
by a power spectrum index n. As mentioned above, both 
indexes appear to be similar for supersonic turbulence, with 
n' ~ n ~ 3.8 (Kritsuk et al. 2007). Although the mass-size 
relation of (evolved) prestellar cores cannot be directly ap- 
plied to the one characteristic of the initial mass reservoirs, 
it is interesting, however, to note that a value n — 3.8 in 
Eq. ([3]) yields a relation for the bound structures supported 
by turbulence M oc R 18 (see Eq.(43) and below in HC08 for 
the complete relation), consistent with the aforementioned 
observational determination from Herschel. 

Interestingly, using the aforementioned observationally 
inferred value D sa 2.5-2.7 yields a = 3.5-3.7, according to 



1 Assuming a Gaussian (lognormal) density probability in- 
tead of a 
of HC08). 



stead of a uniform one yields in that case 4^ oc M 3 (Eq.(37) 



Eq.©, in the case of purely thermal support, which is still 
significantly steeper than the Salpeter value, a — 2.25-2.35 
and 2.5-2.6 in the Burgers and Kolmogorov limits, respec- 
tively, and a = 2.39-2.5 when using the aforementioned 
value for the supersonic turbulence index. 

4. Discussion and conclusion 

The complete calculations, as done in HC08 and HC09, 
yield of course more accurate determinations of the slope 
of the CMF, both in the thermal and nonthermal case. It is 
important to stress that, in the present simplistic approach, 
one assumes a space distribution given by Eq. ([5]), which 
leads to the (D-3) correction in Eq. ©. In the rigorous 
(Press-Schechter based) statistical approach developed in 
HC08 and HC09, the dimension D, both for the (unbound) 
clumps and the (bound) cores is naturally encapsulated in 
the (density and velocity) scaling properties of turbulence 
entering the HC theory. There is no need to explicitly as- 
sume a fractal index, so there is no explicit dependence of 
the Salpeter slope upon such an index. On the contrary, the 
size distribution naturally arises from the calculations with 
(as easily inferred from eqn.(42) of HC08 and the mass-size 
relation (j4]) - detailed in HC08 - derived from the viral con- 
dition): dN/dR oc R-^ w i?~ 3 4 , i.e. D « 2.4, consistent 
with the aforementioned Herschel determinations. In con- 
trast, a pure thermal support (eqn. ([3"| and n — 3) in the 

complete HC theory yields dN/dR oc R 2n - 4 oc R , i.e. 
D = 2. The sometimes advocated dependence of the IMF 
upon the spatial dimension D, and the fractional value of 
the latter, thus simply reflect the (stabilizing) impact of the 
turbulent velocity field on large-scale bound structures. 

The simple dimensional relations derived in this short 
note thus enable us to grasp the crucial impact of turbu- 
lent support, combined with gravity through the virial con- 
dition, on the shape of the IMF. Indeed, based on simple 
scaling arguments, these relations demonstrate that non- 
thermal support is mandatory to yield the correct high- 
mass power-law part of the IMF. Pure thermal support 
would produce a much steeper IMF, leading to a drastic 
lack of massive ( > 1 M@) cores/stars. Interestingly, Eq.© 
yields da/dn < 0, as in HC, and thus predicts that the IMF 
should steepen with a decreasing velocity power spectrum 
index. This could be tested with dedicated numerical sim- 
ulations. It also shows that magnetic field is not needed to 
recover the power-law slope of the IMF. 

As mentioned in Hennebelle & Chabrier (2008), the lim- 
itation of the HC theory is the lack of time dependence. 
Indeed, initial massive mass reservoirs that have been se- 
lected, in our statistical formalism, as being stabilized by 
turbulence might in reality fragment into smaller pieces 
during the cloud's evolution. The impact of such a time 
dependence on the IMF is presently under study. Further 
fragmentation of these primordial cores, however, is un- 
likely to significantly affect the high-mass part of the IMF, 
for two reasons. First of all, given the scale-invariance of 
turbulence, although fragmentation of massive cores dur- 
ing evolution may affect the lower mass part of the IMF, 
by producing new small cores, it is unlikely to drastically 
affect the high-mass part, as turbulence- induced fragmen- 
tation processes will be similar for all masses. Statistical 
explorations of core fragmentation indeed show that the 
high-mass part of the IMF is still barely affected (Swift & 
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Williams 2008). Second of all, observations of massive cores 
suggest that fragmentation is rather limited, most of the 
mass of the core ending up in one or just a few smaller cores 
(Bontemps et al. 2010; Longmore ct al. 2011). Indeed, var- 
ious numerical simulations of collapsing dense cores show 
that radiative feedback and magnetic fields drastically re- 
duce the fragmentation process (e.g. Krumholz et al. 2007, 
Commergon et al. 2010, Hennebelle et al. 2011). 

Together with the results obtained in recent dedicated 
numerical simulations (Schmidt et al. 2010), these simple 
dimensional arguments strongly argue in favor of the high- 
mass part of the initial core mass function, hence of the 
resulting IMF, being mainly determined by the presence of 
hydrodynamical turbulent motions at the early stages of 
star formation. It is essential, however, to clearly under- 
stand what such a "turbulent support" corresponds to in 
the present approach (i.e. in the Hennebelle- Chabrier the- 
ory). As briefly outlined in the Introduction, turbulent sup- 
port here is not to be understood in a static sense, as often 
(most of the time incorrectly) in the literature, by suppos- 
ing that turbulence will remain present in the core and will 
act permanently as an extra source of pressure. Indeed, tur- 
bulence within an overdense region will decay within about 
a crossing time. The mass-scale relation stated by eq.(U]) 
will thus eventually evolve into a thermal relation, as given 
by eq.©, as turbulence dissipates during the collapse of 
the core. As mentioned earlier, however, the aforederived 
relations apply to the initial mass reservoirs set up by 
the density fluctuations within the cloud/clump. Turbulent 
support in that sense means that the random motions of 
the flow provide enough kinetic energy at the early stages 
of structure formation for massive reservoirs to circumvent 
the action of gravity, thus stabilizing these otherwise col- 
lapsing structures. In other words, overdense regions with 
lower masses than the turbulent Jeans mass (see Eqs.(25- 
26) of HC08 and eq.(U]) above) will be dispersed back into 
the flow under the action of turbulence within a turbulent 
crossing time, before they become dominated by gravity. 
Only in regions exceeding the turbulent Jeans mass will 
gravity take over, yielding the collapse of the corresponding 
mass reservoirs, providing the initial seeds for the massive 
gravitationally bound prestellar cores. In that sense, the 
concept of a turbulent Jeans mass should be seen as a sta- 
tistical criterion for selecting the initial pieces of gas that 
will be dominated by gravity, and will thus collapse, from 
the ones that will be dispersed, within a turbulent crossing 
time. 
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Appendix 

A point is worth mentioning concerning the slope of the 
CMF/IMF. Let us reason out in 3D, for sake of simplicity. 
In 3D, the case D — 3 corresponds in the k-space to a 
structured space such that N(k) cx k 3 , i.e. to a space filling 
distribution of fluctuations (maximum compactness) . Since 
the mass associated with the fluctuation of scale R is M = 
pR 3 oc fc~ 3 , one immediately gets 

dN » r 2 

dM* M (?) 

The mass spectrum J7]) thus corresponds to the most natu- 
ral mass distribution, set up by fluctuations obeying purely 
geometrical considerations, i.e. a white noise spectrum of 
fluctuations. It is thus not surprising that mass spectra in 
the interstellar medium exhibit a power-law behaviour with 
an exponent close to 2, namely ~ 1.6 for unbound struc- 
tures such as CO clumps, ~ 2.0 for stellar clusters, and 
~ 2.35 for cores/stars. Therefore, the aim of any theory of 
the IMF is to explain the departure from the naturally ex- 
pected a = 2 slope, which arises from physical, not purely 
geometrical, effects. 
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